Cooper pair sizes in superfluid nuclei in a simplified model 
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Cooper pair sizes are evaluated in a simple harmonic oscillator model reproducing the values of 
sophisticated HFB calculations. Underlying reasons for the very small sizes of 2.0-2.5 fm of Cooper 
pairs in the surface of nuclei are analysed. It is shown that the confining properties of the nuclear 
volume is the dominating effect. It is argued that for Cooper pair sizes the LDA idea is particularly 
inadapted. 



Recent studies have revealed surprisingly small exten- 
sions of Cooper pairs in the surface of superfluid nuclei 
PHD]. Such features are potentially very important in 
pair transfers in nuclear reactions [J]. Though the rea- 
son for the small sizes has been identified in our preceding 
paper [j| to be due to the finite size of nuclei, it is nev- 
ertheless instructive to further elaborate on the underly- 
ing reasons of this behavior. We, therefore, develope a 
simplified model which, however, will keep all the essen- 
tial ingredients for the comprehension of the effect. The 
model consists of a spherical harmonic oscillator (HO) 
potential (without spin-orbit) for the mean field together 
with a realistic treatment of pairing using the Gogny D1S 
force [f| . We will see that such a model quite accurately 
reproduces the results for the so-called coherence length 
(CL), i.e. the size of Cooper pairs, of much more sophis- 
ticated selfconsistent HFB calculations 

The questions we will try to answer are the following: 

i) What is the reason for the existence of such very 
small sized Cooper pairs with extensions 2.0-2.5 fm in 
the surface of nuclei, about a factor 2-3 times smaller 
than the smallest size in nuclear matter at low densities? 
Those values are also much smaller than the ones esti- 
mated from the common believe that Cooper pair sizes 
in nuclei are of about the nucleus' dimension [7], what 
is based on pairing in nuclei being in the weak coupling 
regime. Since in weak coupling CL> b, with b the oscil- 
lator length to be used below, the fact that for a nucleus 
with, e.g. nucleon number ^4=120, b ~ 2.2 fm ~ CL m i n ., 
does it lead to the conclusion that in the surface pairing 
is close to strong coupling? Those small sizes also are 
of similar magnitude as the one of the deuteron, that is 
a bound state. Does it mean that the neutron Cooper 
pairs are locally also eventually in a bound state? Actu- 
ally, this might not be completely surprising, since two 
neutrons are almost bound even in free space and pair- 
ing could help to make them truely bound. The question 



then is whether the small size of the CL's is due to partic- 
ularly strong pairing in the nuclear surface (local strong 
coupling) or whether it is essentially due to the confining 
constraints from the nuclear volume? It will be shown 
that the small sizes are dominantly due to the latter ef- 
fect. 

ii) The minimum of the CL, £(i£), in LDA is about 
at the same density as the one in the quantal case []|. 
Then, is the qualitative resemblence of £,(R) calculated 
from nuclear matter in LDA and the quantal £(i?) a for- 
tuitous coincidence, or is that a manifestation of similar 
pairing correlations in both cases? We will see that the 
quantal behavior of £(-/?) in finite nuclei is very similar 
for nominal and almost vanishing pairing. In the latter 
case one should not talk about coherence length but sim- 
ply of the rms distance of uncorrelated pairs coupled to 
angular momentum L=0 which is entirely determined by 
the single particle mean field wave functions. 

We begin our considerations with the density matrix 
corresponding to one major shell of a spherical HO po- 
tential V{R) = foj 2 R 2 with = AlA- 1 ? 3 MeV 

Pjv = /J \rdm >< nlm\, (1) 

nlm 

where the prime on the sum indicates that it only runs 
over all the states \nlm > contained in the major shell 
N. 

We start out transforming this density matrix into 
Wigner (W) space. W-space, or phase space, is useful for 
certain aspects and furthermore it has a well known ana- 
lytic form for the case of a HO potential where it only de- 
pends on the classical Hamiltonian H c i, = p 2 /2m + V(R). 
The corresponding W-distribution is given by [8J 

Pn\w = f N (H cL ) = 8(-l) N e- 2J ^L$ (^f\ , (2) 
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where L^\x) are the generalized Laguerre polynomials. 

We are now ready to present our simplified pairing 
model. We shall write the W-transform Q of the anoma- 
lous density matrix k(t, r') = 

< BCS\a + (r)a + (r')\BCS > as ( spin singlet wave func- 
tion is suppressed) 



k(R, p) 



y^^ N f N (H cl . 

N 



(3) 



with kn = ujvfAr and un, vn the usual BCS amplitudes. 
Please note that the degeneracy factors are missing in 
Eq. ([3]). This stems from the fact that expression @ 
is not normalised to unity but to the degeneracy of the 
shell N. 

The gap parameters A^v can be obtained from the so- 
lution of a gap e quat ion with matrix elements averaged 
over major shells 
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where E N = (3/2 + N)hu, D N = (N + 1)(N + 2)/2 is 
the degeneracy factor of major shell N, and Vn,n' is the 
shell-averaged pairing matrix element To obtain Vjv,iv we 
start from the state-dependent pairing matrix element Q 



< v)\v\$(i/,i/) >= 

< v, D\v\i/, v' > — < v, v\v\v' , v' >, 



(5) 



where the two-particle states \v, v > are product states 
\v > and \v >. The states \v > are represented by 
single-particle wave-functions (j) v {r, a) = 4'nim,{'i f )'4 ) a and 
the corresponding time reversal states by <f>p(r',a) — 
(— l) 1//2_CT 0* im (r')^_ cr . Averaging over the energy shells 
En and E' N , it is easy to show that in phase space the 
shell-averaged pairing matrix elements read (see [lOj for 
more details.). 
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(6) 



with v v (p) = T)v(p) and v(p) being the Fourier transform 
of the Gogny D1S interaction in the 1 So pairing channel 
Q. The factor r\ serves to adjust the pairing intensity by 
hand. 

In Fig. [I] we give the gap at the Fermi energy Af as a 
function of A. We take rj = 0.85 to compensate for the 
fact that we use the bare mass, m* — m what usually 
overestimates pairing. We see that the typical arch struc- 
ture is recovered. Without any averaging, the gap values 
would depend on the individual single particle quantum 
numbers n, I and Fig. 1 would show an additional fine 
structure. In the present case an averaging over the in- 
dividual substates of one major shell has been performed 
not wyping out, however, the essential quantum features. 



FIG. 1: Pairing gap at the Fermi energy computed using the 
Gogny D1S force as a function of the number of nucleons A 
in an isotropic HO potential. 



We now proceed to the calculation of the CL. Inter- 
preting the anomalous density as the wave function of a 
Cooper pair ( we are aware of the fact that this point of 
view has been debated recently [Hj]), the local rms value 
of a pair is given by [l[ 



! N{R) _ Jd 3 ss 2 K 2 (R,s) 
~D{R~) ~ V / d 3 sK 2 ((R,s) 



f jgfr\d K (H cl .)/d(p/h)\< 



(7) 



Here 2R = r + r' and s = r r' and k(R, s) is the Fourier 
transform of re(R, p) of (3). 

Using Eqs. ((2|) and (j3|), denominator and numerator 
under the square root in Eq.Q can be obtained explicitly 
in the case of the HO potential: 
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where a = 1/b — \Jmuj/h is the inverse HO length, K 
and J are the principal HO quantum numbers of the 
shells, kk and kj the corresponding BCS amplitudes of 
the pairing tensor. 

In the upper panel of Fig. 2 we show for different 
values of n. It is seen that only depends very 

weakly on the pairing strength for 77 < 1, this happens 
for instance around the minimum and the similarity 
with the results of the realistic calculations presented 
in [l[ and displayed again in the lower panel of Fig. 2, 
is striking. In particular our model reproduces the very 
small value of in the nuclear surface of about 2fm. 
For r) > 1, the CL starts to move to lower values in the 
interior. However, the minimum again only is very little 
affected. 

In our model it is now rather staightforward to under- 
stand where this striking independence of £,(R) on the 
intensity of pairing comes from. From Eq. ([5]), we can 
realize that the features of fN F (H c i.), where Np corre- 
sponds to the major shell at the Fermi energy, have a 

width of order ~ tuo (one may check this explicitly for 

( 2) 

some low order Lr N ' polynomials). Since in the case of 
nuclei A^v << Tih), the kn are essentially only active at 
the Fermi level and we approximately have from (3) that 
k(R, p) is proportional to fN F {H c i.)- In the limit n — » 0, 
we have the equality (strictly speaking one should in this 
limit change the name and not call it coherence length, 
since there is no coherence any longer; however, for con- 
venience, we will not change the letter £ nor the name) 
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(10) 



where k Hcl {R) = ^(H cL -V(R)) and p Np (R,s) is 

the Fourier transform of fN F (Hci.) with respect to mo- 
mentum p, that is the density matrix corresponding to 
the Fermi level N — Nf- The latter can be obtained 
from (2) as 



Pn(R, s) 



k 1= n 



r 3/2 
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2 Jl 
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K 1= 



With a rescaling of the relative coordinate s — > 2s, we see 
the well known fact, see e.g. [5|, that the density matrix 
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FIG. 2: (Color online) Coherence length for different 
strengths of the pairing force as a function of the radial dis- 
tance R for a symmetric nucleus with A=120. The dashed 
line correspond to the no pairing limit (top). HFB coherence 
length for the nucleus 120 Sn computed for several strengths 
of the pairing D1S Gogny force (bottom). 



for even/odd N is completely symmetric/ antisymmetric 
with respect to an interchange of relative and c.o.m. co- 
ordinates s and R. From Eq. (fTTJ|) it can be seen that the 
dependence of £(R) on A has dropped out completely 
This stems from the fact that in our H.O. model with its 
degenerate shells, in the limit 77 — > 0, the chemical po- 
tential becomes locked exactly at the Fermi level, i.e. at 
the shell Nf- In general, this is not the case in a Woods- 
Saxon potential where it can happen that the chemical 
potential becomes situated in between two subshells. In 
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the upper panel of Fig. [2 we also show the limiting 
value of the coherence length (broken line) when A — > 0. 
It is clear that this asymptotic form of the CL is very 
close to the other curves, in particular at the minimum. 
Therefore, in nuclear physics, in what concerns the CL, 
we are always almost in the asymptotic limit of vanish- 
ing pairing. In a sense, the closenenss of the CL to the 
corresponding uncorrelated value is one of the most strik- 
ing manifestations that nuclei are in the weak coupling 
regime of pairing. Of course, this should not make us for- 
get that on other quantities nuclear pairing has a strong 
influence. A particularly pertinent example, discussed 
recently 0, Q, is the strong influence of parity mixing 
on the spatial features of the (non-normalised) pairing 
tensor. 

Let us now try to analyse from where comes this typi- 
cal behavior of the CL, i.e. of It raises from R=0 
up to R=l-2 fm, followed by a longer almost linear de- 
scent, passing through a shallow minimum of 2-2.5 fm, 
levelling off at some slightly higher asymptotic value. Be- 
fore coming to this study, let us mention again that this 
behavior seems to be very robust being found in realistic 
HFB calculations in nuclei, see Fig. 3 in Q, in a slab 
geometry [12j . as well as in the present very simplified 
HO model. 

Let us consider the normalised square of the den- 
sity matrix, as it enters the definition of the CL. From 
Eas.(fT0]l and (JTTJ> we obtain 



\p N (R, s)| s 



7T 



Jd3 S | Pjv (R, S )|2 47T 

( E^(-1) N - K ^_% (2a 2 R 2 )L^ W) 



.(12) 



In the particularly simple case of N=l, L^ 2 \x) = | — 

x and L^ 2 \x) — 1, and consqequently Eq. (fT2"|) . after 
multiplying by s 4 , reads 

MR,s)|V 



Jds\ Pl (R,s)\ 2 S 2 

a 4 {2R 2 -s 2 /2) 2 s 4 
4a 4 i? 4 -6a 2 R 2 + 3.75' 



(13) 



One sees that that (12) has an R dependent node at 
s = 2R, a feature which is important for interpreting the 
characteristic behavior of the CL. After integrating (12) 
over s, one obtains for the CL 
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where x = 2a 2 R 2 . Minimization with respect to R im- 
plies that 
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FIG. 3: Coherence length (fm) for systems containing A=12 
and A=28 nucleons as a function of the distance to the center 
R (fm) 

the roots of which are: x\ = 2.5 — \/2.5 = 0.9189 and 
X2 = 2.5 + \/2.5 = 4.0818, x\ corresponding to the max- 
imum and X2 to the minimum. 

In the case of N = 1, let us take, somewhat arbitrarily, 
the symmetric open shell nucleus A = 12 (in order to fix 
the value of w a A~ x / Z ). From the definition of x one 
obtains for the position R and the values of maximum 
and minimum of the CL 



R m 
and 

Rn 



1.0315/m t(R) max = 4.3476/m, (16) 
2.1740/m £(ii) miB = 2.0629/m. (17) 



The coherence length corresponding to A = 12 (N = 1) 
is displayed in Fig. |3l It is very surprising that even 
for this simplest case of N = 1 the essential features of 
the CL are already born out. For instance the minimal 
value is at about 2 fm, as in all other cases, realistic 
ones included On the other hand for N = 0, no in- 
dependence of the CL exists. The constant value of CL 
for 7V=0 is about 2 fm for, e.g. A=4, i.e. the a particle. 
Therefore, one needs at least to go to P-shell nuclei, i.e. 
N = 1, so that in the interior the pairs can extend beyond 
2fm. However, coming close to the edge of the nucleus, 
the rms. value of the pair gets (approximately) back to 
its value it has in the a particle. 

The coherence length for A = 28 what correspond to 
a mid shell nucleus with N = 2 reads: 



4a; 2 - 20a; + 15 = 0, 



(15) 



e(R) = — x 

or 

16a; 4 - 224a; 3 - 
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(18) 
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where again x = 2a 2 R 2 . The minimization with respect 
to R implies that 

64a; 6 - 1088a; 5 + 7248a; 4 - 27168a; 3 
+61340a; 2 - 73348a; + 30435 = 0, (19) 

which yields the only real roots: x\ — 0.7822 and X2 — 
7.5622, x\ corresponding to the maximum and x-i to the 
minimum. The other four roots of (|19[) are complex. 
From the previous definition of x one obtains: 

Rmax = 1.0960/m £(R) max = 6.5700/m, (20) 

and 

R mm = 3.4079/m ((fl) mm = 2.1838/m. (21) 

The coherence length corresponding to this case is also 
displayed in Fig|3l In order to better understand the 
qualitatively similar behavior of the CL for the A=12 and 
A=28 cases, we show in Fig. [4] expression (12) and the 
corresponding one for A=28 as a function of s for various 
values of R. The area below the curves in Fig. [4] directly 
yields the CL. The striking feature is that the scenario is 
qualitatively much the same in both cases, inspite of the 
fact that for N = 2 there are two nodes instead of one 
[l3| . The analysis shows that the two nodes also move 
proportional to R from inside to outside in a similar way 
as for the N = 1 case. We surmise that the behavior stays 
more or less the same also for higher N values. There are 
two asymptotic regimes where the nodal structure in s 
practically does not influence the integrand in s ( i.e. 
(12)) and which are more or less dominated by a single 
bump structure. This is the case for very small i?-values 
as well as for large R- values, approximately from the min- 
imum point of the CL onwards. In between, the behavior 
switches from one regime to the other. This is where the 
CL shrinks about linearly with R. In order to exhibit the 
linear behavior more clearly, we show in Fig. 5 the CL's 
for A=12, 28, 120, and 8000. We scale in that figure the 
CL and the R coordinate by the radius at the classical 

turning point, R t = \J -^jj, given by the intersection of 

the chemical potential // = 46.933Afey (remember that 
with a HO potential, fi is independent of of the nucleon 
number A) with the HO potential and, thus, represent- 
ing the size of the system. It is seen that the different 
curves almost are superposed averaging around a linear 
descent. Only the beginning and the ends vary. The po- 
sition of the minimum ranges between a little more than 
half of Rt for A=12 to about 90 percent of Rt for A=8000. 
In the interior, close to the origine, the pairs occupy the 
whole nuclear volume while approaching the surface they 
steadily shrink to about 2-2.5 fm due to the close presence 
of the confinement. After the minimum, i.e. more or less 
after the classical turning point (for very small systems 
the latter does not have such a well defined meaning), the 
pair wave function enters the evanescent region and again 
slightly expands before reaching the asymptotic value. It 



also is worth mentioning that the minimum of the CL 
is very slowly increasing with particle number, approx- 
imately as ~ A 1 / 6 . For A = 8000 the minimum value 
is about 4fm. Once u> — > 0, the CL approaches infinity 
everywhere. A generic feature also is that, independent 
of A, starting from the center at R — 0, the pairs first 
slightly expand up to R ~ lfm before becoming smaller 
getting closer to the border of the mean field. For these 
i?-values around the origine, the nodes lie in the region 
which is dominated by the phase space factors s 2 and 
,s 4 in the integrals over s in (13), i.e. well to the left 
of the maximum of the bump created by the function 
s 4 e~ a s I 2 . It can easily be verified from our 'easy' ex- 
ample N = 1, eqs (12) and (13), that for very small 
values of i?, the surface corresponding to the s-integral 
of the denominator decreases faster than the one of the 
numerator. Therefore, the CL increases. However, once 
the node comes into the region where the exponential 
regime takes over, i.e. where the extension of the system 
is felt, the CL starts its regression. These considerations 
may be elaborated in all details for the case N = 1 and 
also further be elucidated in considering as a complement 
to the density matrix pjv(R, s) its Wigner representation 
Eq © . Not to make the present discussion too heavy, we 
refrain from entering these more detailed considerations. 
The case N = 1 is, as seen, already characteristic and 
can be studied straightforwardly. 

We also should mention that even in our averaging over 
major shells orbit mixing within the shell takes place. 
The cross terms give raise to a distructive interference 
still lowering the minimum of the CL by a small but 
definite amount of about 0.5 fm from its non averaged 
values. This can be realised in comparing Fig. 6 where 
the CL, i.e. local in R rms radii from individual HO 
orbits are displayed (for a precise definition, see @) with 
the broken line in the upper panel of Fig. 2. 

Intrashcll averaging, therefore, is present even in the 
limit of very small pairing with gap values of the or- 
der of subshell spacings. In Q the same study is per- 
formed with the self consistent HFB orbits, see Fig. 17 
in that reference. It is seen that in the self consistent case 
the reduction of the minimal value of the CL from intra 
shell mixing is about 30 percent and, therefore, some- 
what stronger than in the present simplified HO model. 
Though not completely negligeable, this interesting be- 
havior is nontheless a minor effect with respect to the 
feature we are discussing in this work, namely a surpris- 
ing reduction of the minimal value of the CL by a factor 
3-4 from a simple weak coupling estimate [7j. 

We, therefore, can say that in cases where in finite 
Fermi systems typical values of gap parameters around 
the Fermi energy are smaller than Hlo or energy differ- 
ences between neighboring major shells in spherical nu- 
clei, the size of Cooper pairs in supcrfluid nuclei, or other 
finite Fermi systems, is essentially determined by the spa- 
tial extension of the single particle wave functions close 
to the Fermi energy. Passing for the sake of the ar- 
gument to a continuum version of (2), i.e. k(R, p) = 
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FIG. 4: (Color online) Normalized square of the density ma- 
trix p(R, s) multiplied by s 4 as a function of the relative coor- 
dinate s. In the upper panel it is displayed for A=12 (TV = 1) 
and in the lower panel for A=28 (N = 2) 



J dEK(E)fE(Hd.), we see that in nuclei where the width 
of k(E) is much smaller than the width of fE(H c i.), the 
CL is dominated by fE{H c i.) on the Fermi surface. Of 
course, a very different situation prevails in the opposite 
regime where A >> fuo. In the extreme case of infinite 
matter or hu> — > 0, we have fE{H c i ) — > S(E — H c i ), and 
the ratio of the values of the widths is inversed! Simple 
scaling arguments show that in the latter case £ ~ 1/A 
which also is reflected in the well known expression given 
by Pippard 



_ 1 h 2 k F 
7r m A 



(22) 



or by an equivalent formula given in Appendix of our 
preceding paper [5]. Therefore, in the infinite matter 
case the dependence on the gap is not at all compensated 



FIG. 5: (Color online) Coherence length (CL) for systems 
containing A— 12, 28, 120 and 8000 nucleons as a function 
of the radial distance R. Notice that both, coherence length 
and radial distance, have been scaled by the classical turning 
point distance R t (see text) 
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FIG. 6: Individual contributions to the CL in the HO model 
for A=120 corresponding to the different individual orbitals 
of the last shell. 



between numerator and denominator in eq (6), whereas 
this is the case in finite nuclei, see eq (9). 
As a consequence, the use of LDA, which is equivalent to 
the infinite matter regime, is not valid to estimate the co- 
herence length in finite nuclei. For other quantities, how- 
ever, as, e.g. the pairing energy, LDA gives a reasonable 
good average 15(. Nevertheless, even in such favorable 
cases, LDA is very much at the limit of its validity, for 



instance in what concerns a detailed description of the 
radius dependence of various pairing quantities. Further 
considerations on this subject will be published elsewhere 
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In conclusion, concerning the extension of Cooper pairs 
in finite superfluid Fermi systems, we have identified two 
regimes: One for fkj >> A where the coherence length 
is practically independent of A and determined by the 
spatial extension of the single particle wave functions. 
Besides in nuclei, such a situation may be found in ul- 
trasmall superconducting metallic grains [13, [3 • In the 
second regime with friu << A, the coherence length is 
approximately inversely proportional to the gap values. 
The latter situation is, besides nuclear matter, e.g. re- 
alised in cold superfluid fermionic atoms in traps where 
typical values of A/(haj) may be of the order of ten or 
even larger [l9j . It would be interesting to study the 
cross over from one regime to the other in more detail. 

Let us finally wrap up the situation of the CL in nuclei. 
We found that a simplest spherical HO model already 
simulates quite faithfully realistic HFB calculations with 
the Gogny force. In what concerns the CL, the situa- 
tion for nuclei is such that there is very little difference 
between rms. values of uncorrelated pairs coupled to 
L = S = calculated locally as a function of the radius 
R and local Cooper pair sizes calculated with the nomi- 
nal pairing interaction. Therefore, the small Cooper pair 
size of 2.0-2.5 fm in the surface of nuclei is practically 
entirely a finite size effect and has not much to do with 
existing enhancement of pairing in the nuclear surface. A 
very characteristic and generic pattern has emerged. In 
the lightest nuclei, like, e.g., the a particle, their size is so 
small that the extension of a pair cannot reach more than 
2 fm. Going to P shell nuclei, in the interior the pairs can 
already somewhat extend but approaching the border of 
the mean field they shrink until they again reach a value 
of around 2 fm due to the resricted space around the sur- 
face. In the interior the pairs grow approximately with 
the size of the nucleus, see Fig. 5 but towards the sur- 
face they always regress to their very small value. In the 
evanescent region, the pair sizes become slightly larger 
than their minimum value in the surface region but this 
increase is very moderate. A characteristic feature also 
is that the pairs only feel the finite size from R=l fm 
onwards. Before, they slightly expand up to R=l fm, 
independent of the mass number of nuclei and of parity 



of the shell. This scenario of a first slight increase, fol- 
lowed by a longer linear descent, before going through a 
shallow minimum at 2.0 - 2.5 fm, levelling off in a slightly 
increased asymptotic value is practically a generic feature 
of local pair sizes in nuclei. It is seen in our schematic 
model, but also in realistic calculations, see Fig. 3 of 
though in the latter case some scatter exists, probably 
due to more pronounced shell effects. This characteris- 
tic pattern of local Cooper pair sizes, practically inde- 
pendent of the strength of the pairing force as long as 
it stays below the nominal value, is one of the clearest 
theoretical manifestations that nuclei are in a weak cou- 
pling regime characterised by gap values A << hu>. In 
the opposite limit A >> hu, as prevails in infinite mat- 
ter, but in the regime where A < < /i that is still in weak 
coupling (20j . the coherence length varies inversely pro- 
portional to the gap, a fact which is well known. The fact 
that Cooper pair sizes are largely dominated by geome- 
try should, however, not make us forget that for other 
quantities nuclear superfluidity has an enormous impact. 
To say it again, a particularly striking example, besides 
others more standard ones, is the effect of parity mixing 
on the spatial behavior of the non local (unnormalised) 
pairing tensor, as revealed recently 0,0- Indeed, this 
not normalised pairing tensor k(R, s) becomes very much 
localised in s around the i?-axis whereas the parity pro- 
jected k{R, s) is completely delocalised 0, Q However, 
due to normalisation in the coherence length this feature 
is cancelled out. 
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